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Abstract 

We calculate the shape dependence of entanglement entropy in (5 + 1) -dimensional con- 
formal field theory in terms of the extrinsic curvature of the entangling surface, the opening 
angles of possible conical singularities, and the conformal anomaly coefficients, which are 
required to obey a single constraint. An important special case of this result is given by 
the interacting (2, 0) theory describing a large number of coincident M5-branes. To de- 
rive the more general result we rely crucially on the holographic prescription for calculating 
entanglement entropy using Lovelock gravity. We test the conjecture by relating the en- 
tanglement entropy of the free massless (1,0) hypermultiplet in (5 + l)-dimensions to the 
entanglement entropy of the free massive chiral multiplet in (2 + l)-dimensions, which we 
calculate numerically using lattice techniques. We also present a numerical calculation of the 
(2 + l)-dimensional renormalized entanglement entropy for the free massive Dirac fermion, 
which is shown to be consistent with the F-theorem. 



December 12, 2012 



Contents 

1 Introduction Q] 

2 The shape dependence of EE in (5 + l)-dimensional CFT [4] 

2.1 EE with non- vanishing extrinsic curvature [5] 

2.1.1 EE in the interacting M = (2, 0) theory 

2.1.2 EE using Lovelock gravity [12] 

2.2 EE with conical singularities [13] 

3 EE in (2 + l)-dimensional free massive theories [T6l 

3.1 Numerical computation of the renormalized EE [TS] 

4 Discussion 22 



1 Introduction 

The ground state entanglement entropy (EE) in quantum field theory measures the quantum 
correlations between two subspaces separated by a surface S, called the entangling surface 
(see, for example, [lfj6]). More precisely, if p = |0)(0| is the density matrix of the ground 
state, then we are to construct the reduced density matrix pr obtained by tracing over 
the degrees of freedom within one of the subspaces. The entanglement entropy is defined 
as the Von Neumann entropy of the reduced density matrix: S = — tr(pR log pr). The 
entanglement entropy is dominated by short-range correlations across the entangling surface. 
If e is the short-distance cut-off of the (d + l)-dimensional field theory, these correlations give 
a contribution to the entanglement entropy proportional to the area of S divided by e^ 1 . 

The subleading terms in the entanglement entropy contain useful, cut-off independent 
information about the field theory. One application of these terms has been to characterize 
the topological order of the gapped ground state of many-body systems (7j[8). Another 
application was to propose the F-theorem in (2 + l)-dimensional QFT (9 12 , which states 



that the negative of the finite part of the entanglement entropy decreases under RG flow 
from a UV fixed point to an IR fixed point. In even dimensional CFT it is expected that the 
entanglement entropy across a smooth entangling surface contains a loge divergence, with 
the coefficient of the term related to the conformal anomaly. In odd dimensions there is no 
conformal anomaly, and as a result when the entangling surface E is a smooth submanifold 
the entanglement entropy is not expected to have a loge divergence. One way to see how the 
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loge divergence arises is to consider the replica trick for calculating entanglement entropy. 
In this approach, the entanglement entropy across a surface E is related to the partition 
function of the theory on a copy of the background spacetime which is conically singular 
along the entangling surface. 

The form of the loge divergence in the entanglement entropy is known in 2 and 4- 
spacetime dimensions. In (1 + 1) -dimensional CFT the entanglement entropy across an 
interval of length L takes the simple form Sl = — (c/3) loge + 0(e°), where c is the central 
charge. In (3 + l)-dimensional CFT there are two Weyl anomaly coefficient^] a and c. With 
the entangling surface E embedded in IR 1 ' 3 , the CFT entanglement entropy is given by [13] : 



* = "7? + 1 180 / E * + 2407 X <" - 2* ' j iog£ + o(f ' • ( ") 

where i?2 is the Euler density on E normalized such that j s2 E 2 = 2, and k % ab is the extrinsic 
curvature, with a,b, ... local indices on E and z = 0, 1 labeling the unit normal vectors 
ril, one of which is timelike and the other spacelike. We use the standard notation k 2 = 
^ab^cd^^^ anc j ^ r ^,2 = ^ad^bc^.fci^y^ with the induced metric on S and 77^ the 

2-dimensional Minkowski metric. The functionals of extrinsic curvature which multiply a 
and c are invariant under diffeomorphisms on E and under conformal transformations. 

One of our main results is to find the form of the log e divergent term in the entanglement 
entropy in (5 + l)-dimensional CFT in flat Minkowski space when the smooth entangling 
surface E is embedded in a constant-time hyperplane (k® b = 0). There are 4 Weyl anomaly 
coefficients (A, Bi, B 2 , B 3 ) in (5 + l)-dimensional CFT. Using the 2 and 4-dimensional cases 
as motivation, we hypothesize that the loge divergence in the entanglement entropy will 
schematically be of the form Ss\i og = (B\ fi(k)+B 2 f 2 {k)+B 3 f 3 (k)+A fi{k)) log e, where the 
4 functionals fi(k) are invariant under diffeomorphisms on E and conformal transformations 
which keep k® b = 0. We find a basis of functionals for the fi(k). In the case where the anomaly 
coefficients are related by B 3 = [B 2 — (5i/2)]/3, we can do much better and precisely calculate 
the coefficient of log e in terms of the anomaly coefficients and the extrinsic curvature. A 
special case of this result is found in the interacting Af = (2, 0) theory describing a large 



number N of coincident M5-branes (see, for example, 14-18]), for which all of the anomaly 
coefficients are related in a simple way. 

In the case of the interacting M = (2, 0) theory we use the fact that the theory is dual 
to the low-energy limit of M-theory on AdSj x S4 and utilize the holographic procedure for 

: We normalize a and c so that they both equal one for the real scalar. 
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calculating entanglement entropy [19j-[22] . For the more general case where the only relation 
between the central charges is B% = [B2 — (I?i/2)]/3, we proceed by using the prescription 
in 



23 , 24 for calculating entanglement entropy holographically with Lovelock gravity 25 



When E is allowed to have conical singularities we find new divergences in the entangle- 
ment entropy ranging in severity from 1/e 3 to log 2 e depending on the codimension p of the 
singular locus, with p — 1, . . . , 4. The most singular behavior comes from entangling surfaces 
for which p = 1. The log 2 e divergence occurs when E has a conical singularity located at a 
point within E (p = 4), and in this case we find that the coefficient of the log 2 e term depends 
simply on elementary functions of the opening angle of the cone. Similar behavior has been 



found in (1 + l)-dimensional CFT 26-28 , where E can have cusp (p = 1) singularities, 



and in (3 + l)-dimensional CFT 29 , where E can have wedge (p = 1) and cone (p 
singularities. 



A novel check of Solodukhin's formula in (1.1) was given in 30 , where the coefficient of 
the l/(mR) term in the IR expansion of the entanglement entropy across a circle of radius R 
for massive free scalar and Dirac fields in (2 + l)-dimensions was shown to be related to the 
coefficient of the log e term in the entanglement entropy of the corresponding massless theory 
in (3 + l)-dimensions. The entangling surface in (3+l)-dimensions is taken to be E = S 1 x S 1 , 
where the first circle has a large radius L and the second has a radius R <C L. The author 
then numerically calculated the coefficient of the l/(mR) term in (2 + l)-dimensions using 
the methods developed in (4 31,32 . Perfect agreement was found between the numerical 
calculation and the calculation based on the relation to (3 + l)-dimensional entanglement 



entropy. This provided a powerful check of (1.1). 

We perform an analogous check on our formula for the coefficient of the log e term in the 
(5 + l)-dimensional CFT entanglement entropy. In this case the l/(mR) 3 term in the IR 
expansion of the massive (2 + l)-dimensional free field entanglement entropy across a circle 
of radius R is related to the entanglement entropy of the corresponding massless theory in 



(5 + l)-dimensions 29 . The entangling surface in (5 + l)-dimensions is taken to be T 3 x S 1 , 
where the 3-torus has a large volume compared to R 3 , and R is the radius of the circle. The 
constraint B 3 = [B 2 — (-Bi/2)]/3 in our (5 + l)-dimensional formula means that we must 
consider the free massive chiral multiplet in (2 + l)-dimensions, which consists of two real 
scalars and a Dirac fermion. We find good agreement between the numerical calculation in 
(2 + l)-dimensions and the prediction of our (5 + l)-dimensional formula. 

Additionally, we numerically calculate the renormalized entanglement entropy [33] for 
the massive Dirac fermion in (2 + l)-dimensions. The renormalized entanglement entropy is 



3 



defined in terms of the entanglement entropy S(R) across the circle of radius R by F(R) 
—S(R) + RS'(R). At conformal fixed points it was shown in 



11 



,[34] that T(R) = F, where 
F is the finite part of the free energy of the CFT conformally coupled to the round S 3 . 
For a field theory that flows from a UV CFT to an IR CFT, it was conjectured (9 12 that 
Fuv > -Fir- Casini and Huerta p5] showed that ^{R) should be a monotonic interpolating 
function along the RG flow between the UV and IR fixed points. We verify that this is the 
case for the massive Dirac fermion. A similar check was performed for the free massive scalar 
in 



33 



While this paper was in its final stages of preparation, we learned of the work of 136] , which 



partially overlaps with our results in section |2.2| on entanglement entropy across singular 
entangling surfaces. 



2 The shape dependence of EE in (5 + l)-dimensional 
CFT 

The origin of the logarithmically divergent term in (5 + l)-dimensional entanglement entropy 
is the conformal anomaly (see, for example, 18 , 37]) 

3 



(2.1) 



n=l 



where = —{2/^/—g)5I/5g^ 1 ' is the stress tensor, the B n and A are the anomaly coeffi- 
cients, Eq is the 6-dimensional Euler density, and the I n are geometric invariants constructed 
from the Weyl tensor C^ vpa . For the moment we are taking the background manifold to be 
a general pseudo-Riemannian manifold Ai. In our conventions 
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(2.2) 



and in Euclidean signature j s6 E§ = 2. 

When there is a rotational symmetry in the space transverse to E, the logarithmically 
divergent part of the entanglement entropy can be determined [9] using the replica trick 



5- 
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(2.3) 
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where g^ u = n % n % v . The rotational symmetry requirement forces us to consider entangling 
surfaces with vanishing extrinsic curvature. Hung et al. [23] have proposed that when the 
restriction of rotational symmetry in the space transverse to E is lifted but one maintains 



k l ab = 0, the result in (2.3) is modified by a term proportional to B 3 . The expression in (2.3) 
is invariant under conformal transformations g^ u — > e~ 2w g /J)U for which n^d^oj = 0. These 
are conformal transformations which preserve k l ab = 0. The term that fails to be invariant 
under more general conformal transformations is the one proportional to B%. 



Notice that when the background metric is conformally flat, the expression in (2.3) van- 
ishes. In this case the only contributions to the logarithmically divergent part of the entan- 
glement entropy come from the extrinsic curvature. 

2.1 EE with non- vanishing extrinsic curvature 

We now calculate the dependence of the entanglement entropy on the extrinsic curvature of 
the 4-dimensional entangling surface E. For simplicity we take the background to be flat 
IR 1,5 , and we take vanishing extrinsic curvature in the time-like direction (k° ab = 0). For now 
we assume that the coefficient of the A-anomaly term remains a topological invariant when 
E is allowed to have extrinsic curvature. While this is a common assumption, it will be 



checked holographically in section 2.1.2 With this assumption we have Ss = (Bi fi(k) + 
B2 f2(k) + B 3 fo(k) — A J s E4) loge, where the functional fi(k) of the extrinsic curvature 
should be invariant under local diffeomorphisms on E and conformal transformations which 
keep k° ab = 0. They should also vanish when E is the round S 4 , since the CFT entanglement 
entropy across the round S 4 is equivalent to the Euclidean free energy of the CFT conformally 
coupled to S 6 11 34 . The log e term in the S 6 free energy is proportional to the A- anomaly 
coefficient (see, for example, [38]). 

We need to find 3 appropriate functionals of the extrinsic curvature to span the fi(k). We 
first consider functionals which only depend on the extrinsic curvature and not on derivatives 
of the extrinsic curvature. The most general form of such a functional is 

T = J (ai k 4 + a 2 k 2 tr k 2 + a 3 (tr k 2 ) 2 + a^ktrk 3 + a 5 tr k 4 ) , (2.4) 

where we are using the notation tr k n = k^ a2 k\ 2 a:i ■ ■ ■ k\ n ai , and the coefficients {«i, . . . , a$} 
are to be determined. Under a conformal transformation g^ v — > e~ 2ui g^ Vl the extrinsic 
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curvature transforms non-covariantly, 



kl, = e-»(kl b + ^8,00) , (2.5) 



with -y a b the induced metric on S. The non-covariant factor on the right hand side of (2.5) 
is removed by subtracting off the trace of k\ b . That is, we define the traceless tensor 

k ab = kl b -^fk, (2.6) 

so that under conformal transformations k ab — > e~^k ab . Using the improved tensor k ab we 
may immediately write down the two independent, conformally invariant functionals of the 
extrinsic curvature: 

T 1 = J (^ k2 ) 2 = l ((trA; 2 ) 2 -^ 2 trA; 2 + l^ , 

T 2 = f tr¥= [ (trk 4 -ktrk 3 + -k 2 trk 2 -—k 4 ) . 
k k\ 8 64 J 

The third functional should then involve derivatives of the extrinsic curvature. The most 
general form of such a functional is 

T 3 = J ((V a k) 2 + C!k 4 + c 2 k 2 tik 2 + c 3 (tik 2 ) 2 + dktrk 3 + c 5 trk 4 ) , (2.8) 

where the constants {c\, . . . , C5} are to be determined. We will determine them later in this 
section, but for now we can immediately write down one constraint on the c's. This is the 
constraint which comes from the requirement that T 3 must vanish on the round S A , 

4 3 Cl +4 2 c 2 + 4c 3 + 4c4+ c 5 = 0. (2.9) 

An important observation is that when the background manifold M. is allowed to have a 
nonzero curvature R pV pai the functionals T\ and T 2 remain conformally invariant functionals 
without the need for adding any new terms proportional to the background curvature. This 
cannot be the case for the functional T3. This can be seen quite simply be restricting to 
conformal transformations — > e -2 ^^ such that n^d^ou = 0. In this case k — > e +w k 
acts just like a scalar field in 4-dimensions. It is well known that the action of a scalar field 
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conformally coupled to curvature on a 4- dimensional manifold S is 



(2.10) 



where 1Z is the curvature scalar on S. The curvature scalar 71 is related to the curvature of 
Ai by the contracted Gauss-Codazzi equation, 1Z = R ab a b + k 2 — tik 2 , where R a bcd denotes 
the projection of R^ upa onto S. The action I has a term proportional to the integral over 
R ab ab time k 2 . We can similarly expect T 3 to be modified by terms proportional to the 
background curvature when the background curvature is non-vanishing. 



Recall that (2.3), which gives the entanglement entropy in the case where there is a 



rotational symmetry in the space transverse to S, is only invariant under conformal trans- 
formations for which n^d^bJ = 0. The term that fails to be invariant is proportional to the 
anomaly coefficient B 3 . This leads us to conjecture that in the opposite limit, when k\ h is 
non-zero but R 



0, the functional T3 will enter the entanglement entropy through a 
term proportional to B 3 . This conjecture will be checked holographically later in this section. 

To summarize the above discussions, we have narrowed the form of the log e term in the 
entanglement entropy to 



'S log 



£ ^B)j T X + $>b\ T 2 + ^ B 3 T 3 -A jT 



E, 



loge 



(2.11) 



where the 7 coefficients {fi^\ fif\ fi^} have yet to be determined. When the anomaly 
coefficients satisfy certain relations which allow us to use Einstein or Lovelock gravity in the 
bulk, we can precisely determine these coefficients through holographic calculations. 



2.1.1 EE in the interacting Af = (2,0) theory 

In the case of the interacting Af = (2,0) theory describing a large number iV of coincident 
M5-brane, it is convenient to write ( |2.11[ ) in the form 



[ VlE 4 + 



96tt 2 



loge 



(2.12) 
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where we have used the anomaly coefficients 18 , 37j 



1=^(1) <*"•>. a^d™)^, (213) 

ft = (4^7! (420) (4Ar3) ' Ba = WW. (140) (4/V3) ' 

We may calculate the entanglement entropy holographically by finding the 5-dimensional 
surface £5, which approaches the entangling surface £ at the boundary of AdS?, is extended 
in the rest of the spatial dimensions, and minimizes the area functional 



1 



S* = — m I dVVdetGffi, (2.14) 
where G^ d is the induced metric on £ 5 . We write the AdS 7 metric as 

ds 2 7 = L\ dS dy2 + d ^ dx \ y>0, (2.15) 

y 

with dx^dx" the metric on IR 1,5 , L AdS = 2£ p (7riV) 1/3 , G% ] = 6 G^ / (it 2 L 4 AdS ) , and G$ 1} = 
16tt 7 ^. 

We consider entangling surfaces of the form £ = T 4_n x £„, where each £ n is a deformed 
n-sphere and 1 < n < 4. In the case n = 4 we simply have £ = £ 4 . We will take the volume 
of each torus T 4 ~ n to be L 4 ~ n , with L large compared to other scales. In each case we write 
the metric on IR 1 ' 5 as 

4— n 

dx^dx* = -dt 2 + dx 2 + dr 2 + r 2 (d6 2 + sin 2 Q d£l\_^) , (2.16) 

where dflf 1 _ l is the metric on the unit S n ~ l . We write the spacetime metric this way because 
we will take £„ to wrap the (n — l)-sphere and be given by the curve r = Riff). The function 
Riff) must be smooth, positive and have vanishing derivative at its endpoints. The latter 
condition ensures that our entangling surface is free of conical singularities. The torus T 4 ~ n 
is extended in the Xi directions, which are compactified such that Xj ~ X, + L. We choose 
such entangling surfaces because they are not overly complicated to describe but still provide 



2 Note that our A anomaly coefficient differs from those in (i~8l[37| and also the recent work of [39] because 
we normalize the Euler density E§ so that f se E§ = 2, while those references choose a different normalization. 
For example, in 37 they use the normalization ~ g^jr (4^)3 Js<> ^& = and in 139 they use the normalization 



•3! (il) 3 is 6 E( > ~ 2 - 



s 



non-trivial geometries. 

Before beginning the calculations, note that it is actually sufficient to just consider E = 
T 2 x E2 or E = T 1 x E3 alone. That is, each one of these entangling surfaces is general enough 
to determine the undetermined parameters. We provide an analysis of all four entangling 
surfaces, however, to show that the method is consistent. In addition to the calculations 
presented below, we have also performed numerical calculations where the radial coordinate 
r in E n is taken to be a function of 9 and the n — 1 coordinates in the internal S n ~ l . Those 
calculations were consistent with the general findings below. 

E = T 3 x Ei 



When E = T 3 x E x the functionals T x , T 2 (2.7) and T 3 (2.8) reduce to 



f A m 21 L 3 1 . 

Ti = —r k\ T 2 = — — / k 



16 J tl 64 , Si 



(2.17) 

T 3 = L 3 / [(Vk) 2 + ci ft 4 + c 2 k 2 tr k 2 + c 3 (tr k 2 ) 2 + c 4 k tr k 3 + c 5 tr n A ] , 
where 

R 2 (6) + 2R' 2 (6)-R(6)R"(6) , x 

k = — \; 9 v ; 2.18 

{R 2 {B) + R' 2 {6)f 2 V 7 

is the extrinsic curvature of Ei in M 2 . We take the bulk entangling surface to be specified by 
the function r = r(y, 9). Near the boundary at y — the solution to the equation of motion 
is 

r{y,9) = R{9) y + 0{y ) . (2.19) 



Substituting the solution into (2.14) we calculate the logarithmically divergent term in the 
entanglement entropy, 

s ^=-&iX m2 -^) 10 ^ (2 - 2o) 
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Comparing with (2.8) and (2.12) we find the following constraints among the undetermined 
coefficients, 



(i) 



1 " 12" + 27^ 3(Cl + ° 2) + 5 ° 3 ) ' M = _1 • 



(2.21) 



In writing down the above expression we have made the choices C4 = (8/3) C3 and C5 = —2 C3 
for later convenience. We are allowed to do this because we can always redefine T 3 by adding 
a linear combination of T\ and T2. We will make these choices for c\ and C5 throughout the 
rest of this section. 



S = T 2 x E 5 



When £ = T 2 x £ 2 , the extrinsic curvature tensor on £ 2 becomes 



1 / R 2 + 2rR' 2 - RR" 
V# 2 + R' 2 1 Rsin 2 6{R- cot 6 R') 



(2.22) 



from which we may calculate the functionals Ti, T 2 (2.7) and T3 (2.8). The holographic cal 



culation is analogous to the previous one. On the bulk entangling surface we take r = r(y,9). 
The solution to the equation of motion for the minimal area surface near the boundary is 



r{y, 0) = m + A, ( cot Bit - Ri?Rr2 t R( l F !r H " )) I y> + 0(y« 



8R 2 



R 2 + R' 2 



(2.23) 



which allows us to calculate the logarithmically divergent term in the entanglement entropy 



as a function of R{9). Comparing with (2.8) and (2.12) we deduce the following relationships 
among the undetermined coefficients 



ci = ^ + |, c 2 = -l-2c 3 , ^ = ^ = 0, a*W = -1 
lb 3 



(2.24) 



S = S 1 x S 3 

With S = S 1 x E3 the extrinsic curvature tensor entries fc„h on S3 are the same as those 



in (2.22) for a, b = 1, 2. The new non-zero entry is 

sin 2 6 2 sin 2 9{R— cot fli?) 



33 



VWTR 12 



(2.25) 
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where, using the notation of (2.16), we write = G^+sin 2 8 2 d(p 2 The procedure described 
above for calculating the entanglement entropy holographically may be carried out in an 
analogous fashion. The solution for r(y, 6) near the boundary is 



r(y,e) = R(6) + 



1 / , R(4R 2 + R (3R-R" )) , , 

— - 2cot0# * 5 -]y 2 + 0(y 4 ). 2.26 

8R 2 \ R 2 + R 2 1 



Calculating the logarithmically divergent term in the entanglement entropy and comparing 



with (2.8) and (2.12), we find the same set of equations as in (2.24). 



£ = £ 4 

The final case we will consider is when £ = £4. The extrinsic curvature tensor entries k a b on 
£4 are the same as on £3, which was discussed in the previous paragraph, for a, b = 1, 2, 3. 
The new non-zero entry is 



k 



sin 2 6 3 sin 2 6 2 sin 2 6{R - cot OR') 
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VR 2 + R 2 

where we write d£l\ = d6 2 + sin 2 O^dVt^. In the holographic calculation we find 



r(y,9) = R(9) 



1 / , R(hR 2 + R(AR- R")) . ., 



(2.27) 



(2.28) 



near the boundary at y = 0. Calculating the entanglement entropy we deduce that 



Cl " 16 + 3 ' 



c 2 



l-2c 3 



12 



(3) 



-1. 



(2.29) 



The fact that C3 is not determined in the above analysis reflects the fact that in this theory 



£4 can mix with T3, since B3 and A are related. However, with the constraint in (2.9) we 
are led to take c 3 = —25/16. With this choice T 3 becomes equal to 



(V a kf k* + llk 2 tik 2 -Qitik 2 ) 2 -IQktrk 6 + 12tik A ) , (2.30) 

s \ 16 



and we arrive at our final formula for the shape dependence of entanglement entropy in the 
interacting = (2, 0) theory: 



Sy 



-N 



2 f T 3 

3 7 S ^ 4 + 96^ 



loge. 



(2.31) 
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2.1.2 EE using Lovelock gravity 



In this section we generalize (2.31) to theories with three independent anomaly coefficients. 
The single constraint on the anomaly coefficients is B% = [B2 — (-Bi/2)]/3, which allows us 
to compute the entanglement entropy holographically using Lovelock gravity. In Lovelock 



gravity the entanglement entropy functional in (2.14) is modified to |23 



4G 



(7) 
N 



<Fa\ G 



£ 5 



.(5) 



ind 



1 + -^fooL AdS 



f2 r4 

o Joo-^AdS 



(2.32) 



where 1Z 



5 puVpa 



is the curvature tensor of the induced metric on the holographic entangling 



surface E5. We leave out a Gibbons-Hawking boundary term in (2.32) since it will not be 



involved in our discussion. The parameter is a positive root of the equation 



l = /oo-0-/^U. 



(2.33) 



The gravitational theory has an AdS vacuum with radius L^ds- The central charges of the 



dual CFT are given by 40 



B 1 
B 3 



L\ dS -9 + 26/^ + 51/^ 



9,tiG% ] 



288 



L\ dS -9 + 34/ 00 A + 75/^/i 



8nG 



(V) 
N 



1152 



Bo — — -Di 

2 2 1 



.4 



7T 



3 L\ dS 3-10/^-45/^// 



(2.34) 



8ttG 



(7) 
N 



We may calculate the entanglement entropy in an analogous fashion to the way described 
in the previous subsection for the set of entangling surfaces S = T 4 ~ n x S n given in that 
section. The key insight is that, at least in the examples we use, the solutions to the equation 
of motion do not change near the boundary in the presence of nonzero A and \i. What is 



required, then, is to take the functions r = r(y,9) given in (2.19), (2.23), (2.26), (2.28) for 
n = 1,...,4, respectively, and use them in (2.32) to evaluate the logarithmically divergent 



term in the entanglement entropy. Carrying out this procedure leads uniquely to 



5- 



2 I log 



A I E 4 + 6 7T 

s 



3 Bo 



B l 



J + B 3 T 3 



loge, 



(2.35) 
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where 



J 



5 k 2 



tr k z + (tr k) + 2(fc tr k - tr fc 4 



(2.36) 



and T3 is given in (2.30). We should emphasize that (2.35) is only valid for theories which 
satisfy the constraint B% = [B2 — (i?i/2)]/3. With this constraint there are three independent 
anomaly coefficients. This is the reason why there are only three independent functionals of 



the extrinsic curvature in (2.35), while in a generic theory we would expect four (see (2.11)). 
A simple, relevant example of the above procedure is when £ = T 3 x S 1 , with the S 1 of 



constant radius R. In this case the solution for r(y) in (2.19) near the boundary reduces to 
r(y) = R — (y 2 /8R) + 0(y 4 ). Evaluating the entanglement entropy functional in (2.32) on 
this solution, we find that 



-^AdS-k 71 



2G 



(7) 
N 



dy 



•H (3 — 5/00(2 A + OA,//)) 3(1 - 2/^A - 3/» 



3y 5 

3(3-2/ 00 A + 3/^/i) 
2048 R 3 y 



32 Ry 3 



+ 0(y c 



(2.37) 



which leads to 



9L 3 7r 2 (20 5 2 - 7 
32 R 3 



loge 



(2.38) 



2.2 EE with conical singularities 

It is interesting to ask what new types of singularities the entanglement entropy can have 
when the entangling surface S in (5 + l)-dimensional CFT is allowed to have conical sin- 
gularities. The conically singular entangling surfaces can be categorized into four different 
types depending on the codimension p of the singular locus, with p = 1, ... ,4. We will 
consider conically singular entangling surfaces of the form S = T^ 4- ^ x C^ p \ where T^~ p ^ is 
the (4 — p)-dimensional torus and is the p-dimensional cone of opening angle f2. 

The singular entangling surfaces are taken to lie in the t = hyperplane of flat Minkowski 
space. For each p we write the spacetime metric as ds 2 = —dt 2 + Ylt=i + dz 2 + dr 2 + 
r 2 d£l 2 _ 1 , with —00 < z < 00 and < r < 00. The torus T^ 4_p ) wraps the directions Xi, 
which are compactified such that Xi ~ x\ + L with L large compared to other scales. The 
cone is taken to wrap the unit S^ -1 and be given by the line z = (cot Q,) r. 



Even though (2.35) was derived under the assumption of smooth S, we can still try and 
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use this equation to figure out the form of the new singularities in the entanglement entropy 
in the presence of conical singularities. This method was applied successfully in (3 + 1)- 



dimensional CFT in 29 . Roughy speaking the integrand in (2.35) goes as 1/r while the 



integration measure brings in a factor of drr p ~ l for p > 1. The case p = 1 is an exception 
since the extrinsic curvature of vanishes everywhere except at the singular locus, where 
it diverges. Restricting to p > 1 we can do the integral over r with the short distance cut-off 
r > e, and we find that when p < 4 we get divergences in the entanglement entropy of the 
form / (O) L 4_p '(log e) / e 4_p , with /(fi) some function of the opening angle, plus less severe 
divergences. By arguments directly analogous to those presented in |29|, one can show that 
when p — 1 the most severe, new divergence the entanglement entropy gets is an f(Q)L 3 /e 3 
type divergence. A simple way of heuristically arriving at this result is to regulate the 
singular tip of by a semicircle times a 3-torus. 

We will be more precise in the case p — 4. Writing the metric on the S 3 as dVl\ = 



del + sin 2 #i(d#! + sin 2 2 , 
are 



we find that the nonzero components of the extrinsic curvature 



= r cos Q , ke 2 g 2 = kg ldl sin 2 Q x , k H = kg 2 e 2 sin 2 6 2 ■ 



(2.39) 



Naively using (2.35) then gives the log e divergence 



9 7T 3 cos 2 f2 

S C W llo K 2 e = Tno h 121 B l + 236 B * + ( 7 B l ~ 20 B *) COS ( 2 Q )} W 6 



128 sinfi 



(2.40) 



in theories where B3 = [B2 — (-Bi/2)]/3. A similar log 2 e divergence was found in (3 + 1)- 
dimensional CFT in [29] with conical entangling surfaces. As in that case, the function 



in (2.40) which multiplies log e diverges like 1/Vt as Vt goes to zero, and it approaches zero 



as Q goes to tt/2, where the conical singularity disappears. We find that the first correction 
away from Q = 7r/2 comes from the functional T3 and is thus proportional to B3. 



ScW I lor 2 



g e 



54tt 3 5, [Vt 



71 



on 



71 



log 2 e 



(2.41) 



At small fl (2.40) becomes 



CM I log 2 



27 7^(19 Si -36 -B 2 ) 
Q4SL 



+ 0(Q) log 2 e. 



(2.42) 



We may understand this result in a simple way because at small Q we may approximately 
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decompose the cone C (4) into the union of "cylinder-like" entangling surfaces of the form 
I x S 3 , where I denotes the interval of infinitesimal length dL, and the S 3 has a radius 
R = Lsinf2. Here L denotes the distance between the apex of the cone and the "cylinder- 



like" surface. Using (2.35) we determine that the entanglement entropy across this surface 
is 



27ir 3 (l9B 1 - 36 B 2 )dL 



£|log 



32 sinfi 



L 



log(e/L) . 



(2.43) 



Integrating the above expression over L and taking Q to be small gives exactly the expression 



in (2.42). 



There is a small subtlety in deriving (2.40) which also arises in the (3 + l)-dimensional 
case |29|. Since we cannot take the logarithm of a dimensionfull quantity, we know that the 



loge term in (2.35) must be accompanied by a finite term proportional to logr such that we 



can write loge — > log(e/r). The finite term must be proportional to logr since the radius 



r is the only other quantity with dimensions of length in the problem. In deriving (2.40) 
one is then led to perform the integral f — log(e/r) = — | log 2 e + O(loge). This differs by 
a factor of 1/2 from the naive expectation J e log(e) = — log 2 e + O(loge). 



We now give a more precise, holographic derivation of (2.40) along the lines of that 



presented in 29 . This derivation will lead to the same factor of 1/2 that we motivated 
heuristically in the previous paragraph. We write the AdS-j metric as 



ds 2 _ d V 2 - dj2 + dr2 + r2 ( d62 + sin2 6 dQ2 ) 

y 2 



(2.44) 



with the boundary entangling surface specified by 9 = Q at t = 0. When the bulk is described 
by Einstein gravity (A = jj, — 0) the entanglement entropy functional becomes 



71' 



dr r 



dd 



sin 3 9 
V 5 (r,0) 



yW 2 + ?- 2 (l + (^) 2 ), (2.45) 



where we have taken y = y(r,6). The symmetries of AdS spacetime and the cone 
allow us to take the ansatz y(r,8) = r/g{9). It is convenient to further change variables to 
s = cos 9, with g(s) = g(0), and to consider s as a function of g instead of g(s). Substituting 



the ansatz into (2.45) and making the changes of variables, we find that the entanglement 
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entropy functional becomes the following functional for s(g): 



S C (4) 



7T 



2G 



J^Jdgg^l- s(g) 2 ) y/l - s{gf + g%l + g^(gf . (2.46) 



To find the solution for s(g) near the boundary we need to solve the Euler-Lagrange 
equation at large g. A direct calculation gives 



s(g) = s + 



3 sp 
8g 2 



(2.47) 



with s = cosf2. Substituting this solution into (2.46) allows us to determine the log e 
divergence in the entanglement entropy. This term comes from performing the following 
integrals over r and g: 



dr 
r 



r/e i 

dg - 

90 9 



log 2 e 



+ 0(loge). 



(2.48) 



The lower bound on the r integral and the upper bound on the g integral are determined 
by the UV cut-off e, while go denotes the minimum value of g, whose precise value is not 



needed to calculate the log 2 e divergence. Using (2.48) we find that 

S C (4) | 



9 7T 2 COS 2 tt 



I log e 



16384 G% ] s in ^ 



[31 - cos(2fi)] log 2 e 



(2.49) 



which agrees with (2.40) when A = /i = 0. 



When A and \i are non-zero, we may proceed by performing a calculation analogous to 



the one presented above but using the full entanglement entropy functional in (2.32). A key 



insight is that the solution for s(g) near the boundary is still given by (2.47) at nonzero A 
and n. What remains, then, is to evaluate the more complicated functional on this solution. 



Doing so leads exactly to the result in (2.40). 



3 EE in (2 + l)-dimensional free massive theories 



It is a non-trivial and remarkable fact that we may check (2.35) by comparing with the en- 



tanglement entropy of the (2 + l)-dimensional free massive chiral multiplet. The relationship 
between the entanglement entropy of free (2n + 4)-dimensional CFT, with integer n > 0, 



and that of free massive (2 + l)-dimensional entanglement entropy is outlined in 29 based 
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on the work of (4,30,41 

For (2 + l)-dimensional theories theories with a mass gap of order m, the general structure 
of entanglement entropy across a smooth entangling surface £ is expected to have the form 



S's = a \- f3m£s 

e 



00 

. c -l-2n 

7 + 2^ 



n=0 



rn 



2n+l 



(3.1) 



where the coefficients c^ 1 _ 2n are integrals of functions of the extrinsic curvature and its 



derivatives 29 43 , and 7 is the topological entanglement entropy 7 8. Throughout the 



rest of this section we will consider the free massive theory consisting of no real massive 
scalars and n\/2 massive Dirac fermions, for which 7 = and /3 = —(no + m/2)/12 (see, for 
example, [5 , 44 , 45] ) . 



The coefficients c_ 1 _ 2n are related to the entanglement entropy of the free massless theory 
in (2n+4)-dimensions. If the entropy of the (2n+4)-dimensional theory across the entangling 
surface £2*1+2 = T 2n+1 x £ has the anomaly term 



S- 



(2n+4) 



,(2n+4) 



log 



log(e) , 



(3.2) 



then we can immediately read off the coefficient c^ 1 _ 2n in (3.1) 29 : 



-l-2n 



7r(27r)"(2n-l)!! (2n+4) 

V0l(T 2 ™ +1 ) SS 2n + 2 



(3.3) 



The above formula is slightly modified for fermions. Dirac fermions in (2n + 4) dimensions 
are in a 2 n+2 -dimensional representation, which after dimensional reduction reduces to 2 n+1 



(2 + l)-dimensional Dirac fermions. Thus, the right hand side of (3.3) should be divided 
by 2 n+1 for Dirac fermions. For example, we can use the expression for Sj? +1 ^ in (1.1) to 



write 29 30 



—j^( n o + 3ni/ 2 ) j> dsn 2 . 



(3.4) 



To use our expression for s^? in (2.35), we need the free field anomaly coefficients in 
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(5 + l)-dimensions 



37 



A _ § • 3! / 5 (6) 191 w \ _ -1 (28 (6) 896 (6) 



(3.5) 

(6) A 

12) ' 



where rig is the number of scalars and n[l is the number of Dirac fermions. The minimal 
theory which obeys the constraint B3 = [B 2 — (i?i/2)]/3 is that of the free M = (1, 0) hyper- 
multiplet, consisting of one Weyl fermion and 4 real scalars. Under dimensional reduction to 
(2 + l)-dimensions this theory reduces to that of 2 free chiral multiplets. Taking n = 2ni/ 2 
in (2 + l)-dimensions, we find 

This allows us to write down the IR expansion of the renormalized entanglement entropy for 
a free massive chiral multiplet, 

^)=^(i + l(^W +o(i/ ( mfl)6) )- (37) 

In the following subsection we check the expansion above against a numerical calculation of 
the renormalized entanglement entropy. 



3.1 Numerical computation of the renormalized EE 

We compute the (2 + l)-dimensional entanglement entropy numerically for the free massive 



scalar and Dirac fermion following the prescription in 30 , which is based on the works 



of (4,31, 32 . Using these techniques, the renormalized entanglement entropy was calculated 



numerically in 33 for the free massive scalar. 

For the free massive scalar, one first expands the field into modes of integer angular 
momentum n. The radial direction is discretized into N units. The discrete Hamiltonian 
for the n th angular momentum mode is given by H n = \ t[\ + | (fiiK^ifij, with 7Tj the 



conjugate momentum to 4>i an d i = 1, . . . ,N. The matrix K n has nonzero entries 30 

K? = I + n 2 + m 2 , K = 2 + ~ 2 + m 2 , = = - . 

2 i 1 ^J%{% + 1) 
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From here one constructs the two-point correlators X 



!(iT- 1/2 )y and P 



(7Tj7Tj) = };(K +1 / 2 )ij. To calculate the entanglement entropy across a circle of radius R in 
lattice units, one must then reduce the matrices Xij and Pij to the r x r matrices and 
P[j, which are defined by taking 1 < i, j < r with r = R — |. The entanglement entropy is 
then given by 



with 30 



S a 



tr 



S(R) = S + 2j2S r , 



n=l 



v/Xp^ + log (v^W + ^ 



n n 



(3.9) 



(3.10) 



The calculation for the free massive Dirac fermion proceeds analogously to that for the 
scalar. We write the two-component Dirac spinor as ip = (u,v) T . We again discretize the 
radial direction into an N unit lattice. Expanding the field ip into modes of half-integer 
angular momentum n, we write the discrete Hamiltonian as 



(3.n; 



where the 4 matrices M 11 , M 12 , M 21 , and M 22 can conveniently be combined into the 
2N x 2N matrix = with a? /3 = i ; 2 and fc, / = 1, . . . , X. The nonzero 



entries of M n are 30 



Mt k = (-l) k -m. 



M, 



2fc-l,2fc 



M. 



2fc,2fc-3 



i 

2' 



M 12 : 

n 

jQ-2fc,2fe- 



n + 



— I- 



M. 



2fe-l,2fc-2 



2 ' 
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Ml 



M 



2fc-l,2fc+2 



n + 



2 



2fc,2fc+l 



(3.12) 



The two-point correlator < > f° r the n th angular momentum mode is then given by 
the matrix C n = Q(—M n ), which in practice can be calculated by first diagonalizing — M n , 
taking the step-function of the eigenvalues, and then rotating back. The contribution to the 
entanglement entropy of the n th mode is S n (R) = — tr [(1 — C£) log(l — C r n ) + C' n logC£], 
where C r n denotes the restriction of C n to the first (2r + l)-dimensional block. The entangling 



19 



circle radius R is related to r by R = r + |, and the total entanglement entropy is found by 
summing over the contributions of all the angular momentum modes: S(R) = \ J2 n ^n(R), 
with the factor of 1/2 coming from the well known effect of fermion doubling on the radial 



lattice. We refer the interested reader to 30 for further details of the numerical computation. 



We discretize the radial direction into a lattice consisting of N = 200 points. To minimize 
lattice effects, we only calculate the entanglement entropy for 30 < r < 50 in lattice units. 
For both the scalars and fermions, we calculate the entropy for m = .005 • % in inverse 
lattice units, with % = 0, . . . , 16, which approximately allows us to cover the range .15 < 
(mR) < 4. Finite lattice effects are most important for the modes with small angular 
momentum |46|. Following |33|, for the lowest 10 angular momentum modes we perform the 
numerical calculation on lattices consisting of 200 + 10 -j points, with j = 0, ... ,10, and then 
extrapolate to infinite lattice size. As in [30], we sum explicitly over the first 3000 angular 
momentum modes and take into account higher angular momentum modes by fitting 



1 



logn 



logn 



S n = a 2 ^ + b 2 — =- + a A — + 6 4 — — + a 6 — + b, 



logn 



n' 



7V- 



(3.13) 



and then summing over the rest of angular momentum modes using the approximate fit 
function. 

In figure [T] we plot the numerical results for the renormalized entanglement entropy 




mR 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 



mR 



Figure 1: The renormalized entanglement entropy J 7 for the massive real free scalar (left) 
and massive free Dirac fermion (right). The black curves are the results of the numerical 
lattice computations. The orange curves are the analytic 1 / (mR) IR approximations coming 



from the coefficient c_ 1 in (3.4). The dotted red lines are the m = values F$ and Fp for 



the scalar and fermion, respectively. 

for the massive scalar and Dirac fermion along with the analytic 1 / (mR) approximation 
coming from the coefficient in (3.4). The numerical results are in agreement with the 
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30 



At m 



analytic approximation. The coefficient was checked in a similar way in 
0, the renormalized entanglement entropy for the scalar and fermion should approach the 
values F s = (2 log 2 - 3C(3)/tt 2 )/2 4 « .0638 and F D = (2 log 2 + 3C(3)/tt 2 )/2 3 « .2190 
in 12 , respectively. We are able to confirm this to high numerical accuracy. This check was 



performed in 33 for the massless scalar. 

We may test the l/(mR) 3 term in (3.7) by plotting the function 



(mR) 3 T, 



71 



1 



chiral 



24 mR 



71 

256 



O {l/(mR) 2 ) 



(3.14) 



at large values of (mR), where J^hirai is the renormalized entanglement entropy of the massive 
free chiral multiplet. In figure [2] we plot this function along with the hypothesized asymptotic 



(mR)' 



^Chiral _ 

V 24 mR/ 



0.014 
0.012 
0.010 
0.008 
0.006 
0.004 
0.002 



1.0 1.5 2.0 2.5 



3.0 



mR 



Figure 2: The function (mR) 3 (J^hirai — umR) Pitted at large values of (mR), where J-chirai 
is the renormalized entanglement entropy of the free massive chiral multiplet, which consists 
of two real scalar fields and one Dirac fermion. The black points are the results of the 
numerical lattice computation, which is subject to numerical error at this level of precision. 
The dotted red line is the predicted asymptotic value gfg, which we compute using (3.7). 



value 



256' 



While unfortunately numerical error does become significant at this level of 



precision, the numerical results do seem to be in agreement with the asymptotic value 



256 



predicted in (3.7) 
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4 Discussion 



There are a number of interesting ways in which our work could be generalized. One could 
try and generalize (2.35) to allow for entangling surfaces for which k Q ab 7^ 0. Another obvious 



way (2.35) could be generalized is to allow for theories which don't obey the constraint 
B3 — I (i?2 — \Bi) . One way of doing this might be to use a more general gravitational 
theory. However, it should also be possible to derive (2.35) without resorting to holography. 

Another interesting avenue would be to generalize the (5 + l)-dimensional calculation of 
the entanglement entropy across £ = T 3 x S 1 to (4 + 2n)-dimensions with £ = T l+2n x S 1 
and integer n > 0. The resulting expression evaluated for a multiplet consisting of a free 



massless Dirac fermion and 2 n+2 free real massless scalars would, in conjunction with (3.3), 
give the l/(mR) 2n+1 term in the (2 + l)-dimensional renormalized entanglement entropy 
of the massive free chiral multiplet. Since the massive free chiral multiplet is perhaps the 
simplest massive theory in (2+l)-dimensions, it would be interesting to see if the expansion of 
its renormalized entanglement entropy re-sums in a simple way. Carrying out this procedure 
would require calculating the conformal anomaly coefficients (or some combination thereof) 
of the multiplet consisting of a Dirac fermion and 2 n+2 real scalars in (4 + 2n)-dimensions. 

at 



It would also be interesting to connect these results with recent attempts 39 47 49 



proving the 6-dimensional a-theorem 1381. These attempts all use the methods applied by 



Komargodski and Schwimmer in 50,51 to successfully prove the 4-dimensional a-theorem. 
It might be the case, however, that an entanglement entropy based proof the 6-dimensional 



a-theorem, similar to the proofs given in 2 and 3-dimensions by Casini and Huerta 35 , 52 
is also possible and perhaps more tractable. 
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